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We extend our previous calculation of the quasi-local contribution to the self-force on a scalar 
particle to general (not necessarily geodesic) motion in a general spacetime. In addition to the 
general case and the case of a particle at rest in a stationary spacetime, we consider as examples 
a particle held at rest in Reissner-Nordstrom and Kerr-Newman space-times. This allows us to 
most easily analyse the effect of non-geodesic motion on our previous results and also allows for 
comparison to existing results for Schwarzschild spacetime. 



I. INTRODUCTION 



In previous work [T] (Paper I), we calculated the quasi-local contribution to the self-force on a scalar particle in 
geodesic motion in a general spacetime. We now show how this can be extended to allow for the motion to be 
non-geodesic. A strong motivation for doing so is that it allows our results to be compared against existing work. 
Anderson and Hu [2] have calculated the coordinate expansion of the function V{x,x') appearing in the Hadamard 
form of the Green's function, providing a convincing check on the validity of the relations between covariant and 
coordinate expansions given below. In fact, the comparison allowed a missing factor of 2 in the results of Ref. [2\ to 
be discovered [3j ^ . Anderson and Wiseman 5 built on this work to calculate the quasi-local contribution to the 
scalar self- force for the case of a particle held at rest in the Schw arzschild spacetime. This is useful not only as a 
check on the current work (in particular, with Q = 0, for Eqs. (4.3 1) but also allows us to correct their result. Finally, 
Wiseman [Bl has shown that the total self-force in this case of a static scalar particle in Schwarzschild spacetime is 
zero. This could potentially facilitate a study of the usefulness and accuracy of the matched expansion approach [Tl |5] 
to the calculation of the self- force. 

In Sec. |TT] we extend the general equations of motion calculated in Paper I to allow for non-geodesic motion. We 
then show in Sec. |III| how these expressions simplify significantly if the space-time is assumed to be stationary, and 
that they simplify further under the assumption that the space-time is static. As examples, the cases of a particle 
at rest in Reissner-Nordstrom and Kerr-Newman space-times are considered in Sees. |IV| and |V] respectively. The 
Reissner-Nordstrom result is also used as a check on the existing results of Refs. |S] and [5]. 

Throughout this paper, we use units in which G = c = 1 and adopt the sign conventions of [7]. We denote 
symmetrization of indices using brackets (e.g. (a/3)) and exclude indices from symmetrization by surrounding them 
by vertical bars (e.g. (a|/3|7)). Roman letters are used for free indices and Greek letters for indices summed over all 
space-time dimensions. The Roman letters k, I are used for indices over spatial dimensions only. 



II. THE QUASI-LOCAL SCALAR SELF-FORCE AND EQUATIONS OF MOTION 

In Paper I, we calculated an expression for the scalar self-force on a scalar charge q and mass m travelling on 
a curved background spacetime. However, this expression was only valid provided the particle's path was that of 
a geodesic of the background spacetime. Fortunately, with some care, this result can be extended relatively easily 
to allow for non-geodesic particle motion. As the calculation is largely the same as for geodesic motion, we briefly 
review it here (focusing mostly on the differences caused by the non-geodesicity of the motion) and direct the reader 
to Paper I for more extensive coverage. 

We begin with our expression from Paper I (and previously given without some terms in Refs. W and [S]) for the 
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self-force in terms of a sum of local and non-local parts: 

r = Q (a"^ - a^u"") + ^ (R'^'up + Rp^u^u^u^) + QmL^ " ^ (1 " 6^) + Vim J X/^Gret {x,x') dr' 

(2.1) 

where u° is the 4-velocity of the particle, a is the 4-acceleration, a = §^ is the derivative of the 4-acceleration with 
respect to proper time, Gret{x, x') is the retarded scalar Green's function, mficid is the field mass and ^ is the coupling 
to the background scalar curvature. 

The first terms are all local to the particle's position and can be computed immediately without posing any real 
difficulty. Our task is therefore to elucidate the non-local integral term. To this end, we will work with an expression 
for the self-force which just contains the non-local integral term, with the understanding that the local terms can 
easily be added back in later if necessary: 

/Sl = lim I ^''Gret {x, x') dr' (2.2) 

Note that in the specific cases considered in Paper I, we only had geodesic motion {a°- = 0, a = 0) in Ricci-flat 
spacetimes {Rab = 0, i? = 0) and without field mass (wfioid = 0), so that the local terms were all identically zero. 
However, for non-geodesic motion, the 4-acceleration will always be non-zero so there will always be at least those 
terms to be added back in. 

As in Paper I, we will focus only on the the quasi-local contribution to this integral and leave the remaining portion 
to be computed by other means. The Hadamard form for the retarded Green's function is pUl [TT| : 

Gret {x, x') = (x, x') {U (x, x') 5 (ct (x, x')) - V (x, x') 6 (-cr (x, x'))} (2.3) 

where 0- {x,x') is analogous to the Heaviside step-function (i.e. 1 when x' is in the causal past of x, otherwise), 
S {a {x,x')) is the standard Dirac delta function, U{x,x') and V{x,x') are symmetric bi-scalars having the benefit 
that they are regular for x' x, and a {x, x') is the Synge [51 HH US] world function. Using this form for the Green's 
function, the expression for the quasi-local contribution to the scalar self-force simply becomes 8 : 

fQL^-q^ r V''V{x,x')dT' (2.4) 

Jt-At 

where r — At is a matching point chosen so that x(t) and x'(r') are within a convex normal neighborhood and so 



that the remainder of the integral in Eq. (2.2 1 may be evaluated by other means. 

We now expand V{x, x') in two different ways. First, we express it in the form of a covariant Taylor series expansion, 
i.e. an expansion in increasing powers of the derivative of the Synge world function, cr" = Vcr — 0(cr^/^): 

V (x, x') = ^-|^t^oi...a, (x)a"i (x, x') . . . a"" (x, x') (2.5) 

where explicit expressions for the Wai...ap(a;), up to O (a^^^^ are given in Paper I for massless fields in vacuum 
spacetimes. Paper I also describes how they may be obtained for massive fields in general spacetimes from the results 
Ref. (14|. 

In Paper I, we were able to use this expression for V{x,x') to compute the self-force for geodesic motion. Unfortu- 
nately, things are less straightforward when non- geodesic motion is allowed for. The problem arises as a result of the 
presence of ct" in this expression. As demonstrated in Fig. ([ij, it encodes the proper time, Tg, along a geodesic of the 
background space-time through a'^''aa — —Tg^- However, the integral in Eq. (2.4 1 is along the world line of the particle. 



In the previous case of geodesic motion, this was not a problem as in that case Tg is a natural parameter along the 
world line. For non-geodesic motion, this is no longer the case. To proceed with the calculation using this expansion 
of V{x, x') would require us to first express the geodesic proper time Tg in terms of the integration variable, i.e. the 
particle's proper time t. This is a non-trivial task for general motions in general space-times. An easier resolution 
of this problem comes from expressing V{x, x') in a second form, as a non-covariant Taylor series expansion in the 
coordinate separation of the points, Ax°: 



OO 



V (x, x') = }_J ^^^)„,...„^(x)Ax"^ . . . Ax"- (2.6) 
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FIG. 1: The world-line of a static particle (x' - 
proper time, Tg, and the particle proper time, 
while a° is along the time-like geodesic. 



: constant) in Reissner-Nordstrom space-time is not a geodesic, so the geodesic 
r, are not the same. The integral in Eq. (2.4 1 is along the particle world-line, 



where Ax" = x" — x" and the quantities Va-^,,,ap{x) may be expressed in terms of combinations of the (known) 
quantities Wqi...Qp(x) as follows. 

First, the world function a (x, x') may be expressed in terms of an expansion in powers of the coordinate separation 
of X and x', Ax": 



a = ^g„flAx"Ax'' + A^a^Ax'^Ax'^ Ax'^ + B^a^sAx^Ax'^Ax'^Ax^ + CaB^SeAx°'Ax'^Ax'''Ax^A. 



(2.7) 



where the coe fficients Aap-y, Ba/s-ys, Cap-ySe, ■ ■ ■ are to be determined. In order to determine these coefficients, we note 
that Eq. (2.7 1 implies the coordinate expansion of the derivative of a is: 



CTa =gaaAx" + (ga;3,a + 3A„a/3)Ax"Ax'' + (^c</37,a + 4Saa/3^)Ax"Ax''Ax^-H(B„0^5,a + 5C'aa/375)Ax"Ax''AxT'Ax''- 



Substituting expansions (2.7) and (2.8) into the defining relationship 

2(7 = (r'^da 



(2.8) 



(2.9) 



and equating powers of Ax", we get expressions for each coefficient in terms of the lower order coefficients. The lowest 
terms are given by: 



■^abc 
Bated 

Cabcde 



.9{ab,c) 



X 3 9 

A(abc,d) + ( Q9{ab,\a\gcd),P + ^g{ab,\a\J^\l3\cd) + 2^a(a6^|/3|c 



(2.10a) 
(2.10b) 



^ (^B(^abcd,e) + g°'^ (^'^Aa(abB\i3\cde) + '?>Aa{abAcde),l3 + '^g(ab\,a\B\p\cde) + 2^(afcc|,a|5de),/3^ ^(2.10c) 



This procedure may be easily extended to higher orders using a computer algebra package. 

Next, we obtain a relation between the known coefficients of the covariant expansion, Eq. (2.5 1, and those of the 
coordinate expansion, Eq. (2.6), by first substituting Eq. (2.8) into Eq. (2.5) and then equating the two expansions. 
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In this way we find the following expressions for the Vai...ap{x) in terms of the Va-^...ap{x): 



V 
Va 
Vab 

Vabc 
Vabcd 



'^abcde 



Vab + VaT"^ 

Vabcd + Qv^iabT'^a) + 12^^a/3 (^\^iab^% + 25(,"A 



bed) 



2 ,(ft^>cde) 2^'^{ab.cde) ~^ 2^'^{abcd,e) 



(2.11a) 
(2.11b) 
(2.11c) 

(2.11d) 



(2. lie) 
(2.11f) 



where the F"^ are the ChristofFel symbols of the second kind. 

Although Vabcde could alternatively be given in terms of the Wq,i...c[ (x), we have found that the expression we give 
proves easier to work with. It is found by taking five symmetrized partial derivatives of the equation 



V{x,x') = V{x',x) 



(2.12) 



and then taking taking the coincidence limit x' — > x. It is a special case of the general result that follows from taking 
any number of symmetrized partial derivatives: 



1 f \ 

Va,a2...a, = ^ XI (l) i''^)''^{aia2...a,,a,+^...ap) for p odd 



(2.13) 



k=0 

Applying this identity recursively, we can re-express this with all odd lower order coefficients eliminated: 



Uj p_fci 1 Bp_k+lVia,a,...a,^a, + ,...a,) for p odd 



(2.14) 



k even 



where the Bn are the Bernoulli numbers [15] . Thus, these identities determine all odd coefficients in terms of derivatives 
of lower order even coefficients. 

Now, we simply substitute expansion (2.6 1 into Eq. (2.4 1 and, since V{x,x') is a scalar, take a partial rather than 
covariant derivative to get an easily evaluated expression for the self- force (in this case, it is most natural to work 
with an expression for the self- force in covariant rather than contravariant form): 



J a 



QL 



-q 



V,a - Va - {Va,a " Vaa) Acc" -|- ^ {Val3,a " VafSa) Ax" Ax^ 



{vaPjM - VaUja) Ax^Ax'^Ax'^ + ^ {vap^s,a " Va0jSa) Ax" Ax'' Az^ Ax'' + O (Aa;5) dr' (2.15) 
Finally, we obtain the equations of motion by projecting orthogonal and perpendicular to the particle 4-velocity: 

(2.16) 

-Uu" (2.17) 

Ul 

where 

(2.18) 



dm 
~d^ 



pab ^ gab ^ ^a^b 



is the projection orthogonal to w". 

This expression may be evaluated for a specific particle path in a specific spacetime by writing the coordinate 
separations Ax°' in terms of the particle proper time separation (r — r'). Specific examples of such an evaluation are 
given in the following sections. 
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III. PARTICLE "AT REST" IN A STATIONARY SPACETIME 



The expression for the quasi-local contribution to the self-force given in Eq. (2.151 and the subsequent equations of 
motion take on a significantly simple form if the space-time is assumed to be stationary and if we assume the spatial 
coordinate of the particle to be fixed. Stationarity allows us to introduce coordinates (t, x*) such that the metric 
tensor components, gab are all independent of the time coordinate. In addition, for a particle "at rest", = constant, 
and only the time component of the contravariant 4- velocity, u*, is non-zero. These conditions hold, of course, for the 
case of a particle held at rest in Schwarzschild spacetime which has already received much attention in the literature 

Since the spatial coordinates of the particle are held fixed, the points x and x' are now only separated in the time 
direction. Furthermore, in this case it is straightforward to relate the time coordinate to the proper time along the 
particle's world line. As a result we can rewrite the coordinate separation of the points x and x' in terms of the proper 
time: 



Ax' = At = u\T-T') = 



1 



ir-r'), 



Ax' = 0. 



(3.1) 



Substituting this into Eq. (2.15 1 and performing the straightforward integral of powers of (r — r') gives 

1 



J a 



QL 



-q 
1 

'4! 



(-C.a - Va) At ■ 



^ {Vt,a - Vta) U*At^ + ^ {vtt,a " Vua) (u*)^ At^ 
{Vtu,a - Vttta) (u*)'Ar4 + 1 {vutt,a " Vtttta) (u*)"At5 + O (At^^) 



We can now use Eq. (2.1 If I for Vatcde, along with the analogous equations for Va and Vabc 

1 . 



1 



Vabc 



y,(afcc) 



'^{ab,c} 



(3.2) 

(3.3) 
(3.4) 



to eliminate several terms in this expression. While this did not previously prove particularly beneficial in the general 
case, the fact that partial derivatives with respect to t of these fundamentally geometric objects vanish in a stationary 
spacetime means that many of the extra terms introduced by the substitution will also vanish. Indeed it is easy to 
see from Eq. (2.131 that (in this case of purely time separated points) any term of order 2fc -|- 1 can be related to the 
order 2k term: 



1 



Vtt...ttb 



-Vu...tt,b 



(3.5) 



As a result, Eq. (3.2 1 may be taken to arbitrary order to give: 

A 2fc+l/ t\2k 
,Vtt...tt.aAT ^ (U ) ■ 



J a 



QL 



1 



2 ^ (2fc- 



1) 



1 



2fc/_t\2fc-l 



^(2^)!^^^^^^" 



k 



(3.6) 



To proceed further, we benefit from differentiating between the time and spatial components: 



OO 



1 



2fe/„,t\2fe-l 



■vu...uAr'''{u') 



1^ 



"'^ l2^j2fc + l)!''ii-"'' 



A^2fe/„,t^2fc-l 
Vtt...ttiAT [u ) 



(3.7) 
(3.8) 



From this, it is clear that the t component of the quasi-local self-force only appears at even orders in At while 
the spatial components may appear at both even and odd orders. Note, however, that for a static spacetime - i.e. 
imposing time reversal invariance in addition to stationarity - the second term in Eq. (3.8 1 will vanish since each of 
the Vt...u must be zero in that case (as we have an odd number of t's). As a result, in a static space-time, the spatial 
componen t of th e quasi- local self- force will only appear at odd orders in At. 

By Eq. (2.17) and the fact that only the time component of the contravariant 4- velocity is non-zero, we can now 
rewrite the rate of change of the particle's mass as: 



6 



dm 
d7 



2 1 - A 2fe/ t\2 

^' — 1 2 k 



(3.9) 



Similarly, by Eq. (2.16) we can now write the time and spatial components of the mass times 4-acceleration; 

ma%^ = -^-^h+g''h (3.10a) 
9u 

(3.10b) 



ma 



QL 



where ft and Ji are as given in Eqs. (3.7 1 and (3.8 1 



Again, it is interesting to note the effect of imposing that the spacetime be static. For a static spacetime, the metric 
components odd in t vanish, so our result simplifies to: 



= 



fc=0 



(3.11a) 
(3.11b) 



IV. PARTICLE AT REST IN REISSNER-NORDSTROM SPACETIME 



We now look at a specific example, where we calculate the quasi-local contribution to the scalar self- force on a static 
particle (i.e. held at rest) in Reissner-Nordstrom spacetime. We chose this spacetime as while, through its spherical 
symmetry, it retains much of the simplicity of Schwarzschild spacetime, it also illustrates better the difference non- 
geodesic motion makes. This is because its non- vanishing Ricci tensor means that the self-force appears at lower order 
than was seen in Paper I, where we concentrated on vacuum spacetimes. We therefore have more orders in which to 
see the difference between using the covariant (Eq. (2.5 1) and coordinate (Eq. (2.6 1) expansions for V{x,x'). 

The Reissner-Nordstrom spacetime has line element 



1 



2M 

r 



dr^ 



{d0^ 



1 



2M 
r 



dt^ 



(4.1) 



Since this is an example of a particle at rest in a static spacetime, we may use Eqs. (3.9 1 and (3.111 in order 



to calculate the equations of motion. Computing and substituting in the expressions for the relevant vt...t{x) in 
Reissner-Nordstrom spacetime, along with the time component of the 4- velocity. 



2M 
r 



-1/2 



(4.2) 



we arrive at our result for the quasi-local contribution to the equations of motion of a scalar particle held at rest in 
Reissner-Nordstrom : 



ma^L 



(5Q2 - 8Mr -I- 3r2) Ar^ 



1 



;Y^(1568Q*^-7154MQ^ 



pQ2 (g2 _2Mr 
60rii 

- 10647M2g*r 

- 1236Q''r'' - 171 A/V^ - 1566MQV^ + iQA'Pr^ + 144gVS) Ar^ + ©(Ar^) 



6720ri 

2604QV2 - 5265M3qV3 - 7424,MQ'^r^ + 198M''r'^ + 5147M2qV4 



ma, 

ma:, 



QL 

t 

QL 



dm 



= 
= 
= 



(4.3a) 
(4.3b) 
(4.3c) 
(4.3d) 



(Q2 _ 2Mr + r2) 
2(^8 



A7 



1 



1344r 



^ (196Q^ - 588MgV + 429Af2g2^2 _^ 2Q4.Q^r'^ - 18Af 



- 268Mg^r'' + QM^r'^ + 36gV*) Ar* + 0{At'') 



(4.3e) 
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In the limit Q — > 0, this reduces to the Schwarzschild case, which may be compared to Ref. [S]. However, Anderson 
and Wiseman use an incorrect expression to relate Ai to (r — r') and also use the results of Ref. [2] prior to the 
corrections given in the subsequent errata [3J|3]. Once these two issues are corrected, we find our results are in exact 
agreement. 

An alternative approach in the special case of a static, spherically symmetric space-time is to use the Hadamard- 
WKB approach developed in Ref. [2] to calculate the coordinate expansion of the retarded Green's function, i.e. 



V{x,x')^ J2 irjk{t-t'f {cos J ~iy {r^r')' 

i,j,k=0 



(4.4) 



where the coordinate 7 is the angle on the 2-sphere between x and x' . 

Upon doing so for Reissner-Nordstrom spacetime, we find that the coefficients of the expansion of V{x, x') relevant 
to the current static particle calculation arc"'^: 



VOOO = Vqqi = 



and 



vioi 



V200 = 



V201 



g2 {Q^ - 2Mr + r^)h^ 
20rio 

g2 (Q2 _ 2Mr + r2) (SQ^ - 8Mr + Sr^) 



(Q2 _ 2Mr + r^j 

1344ri6 
+ 91/^3 (r-2M) 
(Q2 _ 2Mr + r^) 



20rii 
2 



196Q^ + 12r(17r - 49M)Q'' + {429M^ - 268rM + SGr^) 



(4.5a) 

(4.5b) 
(4.5c) 

(4.5d) 



1568Q^ + 14Q^r(-511A/ + 186r) + 9MV (22A/2 - 19Mr + 4r 



1344ri7 

+ Q^r^ (10647 Af2 - 7424Afr + 1236r2) + QV^ (-5265Af3 + bUlM'^r - 1566Mr^ + lUr^) 1 (4.5e) 

Additionally, for a particle held at rest in Reissner-Nordstrom, we can immediately relate the coordinate point 
separations Ax" to r — r'. 



Ar = 0, A0^O, A0 = 0, At ^ { 1 - 



2M Q2 



-1/2 



(r-r') 



(4.6) 



Taking the partial derivative of Eq. (4.4), relating the coordinate separations to proper time separations using 



Eqs. (4.6) and performing the integral over r', we arrive at a result which is in exact agreement with Eqs. (4.3) 



V. PARTICLE AT REST IN KERR-NEWMAN SPACETIME 

To illustrate the flexibility of the approach presented in this paper, we will now investigate the case of a particle 
in Kerr-Newman spacetime at rest relative to an observer at spatial infinity. We choose Kerr-Newman for analogous 
reasons to those of Sec. |IV| - the non-vanishing Ricci tensor means that the effects of the non-geodesicity of the motion 
are more apparent. While it is not clear whether the WKB approach of Ref. [2 could be extended to such a spacetime, 
the method presented here is easily adapted to any spacetime, including Kerr-Newman. 

In Boyer-Linquist coordinates, the Kerr-Newman metric is [16 : 



A-a^ + z^ „ , 2(A - r2 - a^){a^ ~ z^) ^^^^ , p\ 2 , 2,^2 



-dtd(j)+'—dr^+p^de^ 



{{a^ + r^f ^ A{a^ - z^)) d<t>^ (5.1) 



^ Note here that our definition of the retarded Green's function differs from that of Ref. '(2[ (including the subsequent errata OUI) by a 
factor of 2 in the term involving v{x,x'). The coefficients given here are consistent with the definition of the retarded Green's function 
given in Eq. 1 2.3 1. 
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where A = — 2Mr + + Q^, = r"^ + a? cos^ 9 and z = a cos 0. Here, Q is the charge, M is the mass and a is 
the angular momentum per unit mass of the black hole. 

For the case of a particle in such a spacetime at rest relative to an observer at spatial infinity, the condition on the 
motion is that the spatial components of the 4-velocity vanish, i.e. — 0. Additionally, the time component of the 
4-velocity may be written as: 



^ 



(5.2) 



Now, since Kerr-Newman is an example of a stationary spacetime, we may use Eqs. (3.9 1 and(3.10l in order to 
calculate the equations of motion: 



dm 

d7 



(ea^ + {q2 + r[-2M + r]}' - 4 {Q^ + r[-2M + r]] + z^ + 6a' {Q' ~ 2Mr + r' - z']) 

2 

- ^^^l^2pii (28(9^° {Tr^ + z^} ~ 4Q** { [343M - U%r\r^ + r[7M + 2Mr]z' + 1312*} + {r"* [3565M2 

-3044Mr + QAQr'^] ~ 2r'' [215 A/^ - 2654Mr + IQlQr'^] + [37M2 + 2304Afr - llOOr^] z^ + 684z^} 

- Q* {[2A'/ - r\r^ [2043^/^ - 1532Mr + 276r2] + [~1124M3 + 9971 AfV - 7516Afr2 + 1440r3] z^ 
+r [166A/3 + 3119Af2r - 3140Afr2 + 864r^] z* + [l37Af2 + 2428A//r - 1008r2] z^ - 156z*} 

+ {r^h2A/ + rf [U7M' - 268Afr + 36r'] - 4[2Af - rjr" [85M^ - 739A/V + 447Afr2 - 63r^] z' 
+ Ar' [55Af* + 371Af V - bibM^r' + 327Mr^ - 72r^\ z^ + 4r [109^/^ + hX^M^r - 427 ATr^ + 72r^\ z^ 
+3 [27M2 - 104A/r + 84r2] z^ - 36zi°} - 9A/2 {[2M - r] + 9r2[-2A/ + rfz' + 9[2Af - r\rz'^ + z^} p'' 
+ 120*^^2 |i54g4 _,. i5M2p2 + 30Q2 y_-jMr + 2p2] } + 6a* {42Q'5 [l3r2 + 92^] - 8Q* [(224M - 89r)r3 
+5(28Af - 9r)rz2 + 44z*] - 4571/^ [2A/r - + 7?\ p* + (l459r* + 846r2z2 + 59z*) 

-20Mr (53r2 - 24z2) igQ (r^ - z^) p*] } + 40^ {7Q® [59r2 + 23z2] - [(2135A/ - 906r)r^ 
+r(623Af + 288r)z2 + 690z*] + Q* [r* (3639Af2 - 299AMr + GOlr^) + (670M2 + 1236A/r - 543r2) z' 
+ (55Af 2 + 2214Afr - 1083r2) z* + Glz*^] + 27Af ^ [r'{^2M + - 3r(-2Af + r)z' + z*] p* 

- Q2 [4Af3^ |-524^4 eir^z^ + 41z*) + M' (-2446r^ + 1251r*z2 + 1776r2z* + 95z^) 

+15A/r (61r* - 118r2z2 + 3z*) p^ - 108 (r* - 3r'z' + z*) p*] } ) Ar* + O(Ar^) 



At' 



(5.3a) 



ma = — - 



g2 (Q2 _ 2Afr) (a2 - z^) 
75600A5/2pi7 



11340Q2p®A (2a2 + Q' - 2Mr + - z^) Ar^ 



{23807r* + 36895r2z2 + 9632z*} - {[119701Af - 49004r]r^ + 8[21670Ar - 4053r]j 



.3 .,2 



-r[32923Af + 60726r]z* + 44146z'^} + 2025A/2 {r''[-2M + - 3r[-2A/ + rjz^ + z*} 

- 3Q* {2Ar2 [32985r6 + 44155r*z2 + 5033r2z* + 775zS] - QMr [8618r* - 3695r2z2 - 7549z*] p^ 

- [9749r* - 20300r2z2 + 2999z*] p*} - 3Q2 [M^ [37174r^ + 47910r5z2 + 5970r3z* + 4450rz6] 

- [40349r^ - 20415r*z2 - 34905r2z* - 2725z6] p^ + 9Mr [l509r* - imQr'^z'' + 295z*] p* 
-1350 [r* - 3r2z2 + z*] p^} + 675a*p* {124Q* + 15Af2p2 + 30g2 [-5A/r + p^] } 

- 675a2p2 {q6 ^^45^2 _^ ^Q3^2j _ 2Q* [(227A/ - 84r)r3 + (143A/ - 30r)rz2 + 54z*] - 15Af2 [2Mr 



^] p* + 2Q2 [m^ (7r2 + z^) (25r2 + llz^) - 2A/r {bQr'^ ~ iiz'^) p' 



15 (r 



)p*]} Ar* + 0(Ar6) 



(5.3b) 



ma 



q' {a' + Q' + r{-2M + r)) 
6720A2pi8 



112QV^A ( 30aV + r {Q' + r[-2M + r] ] {hQ^ + r[~%M + 3r] } + 2 {-2M' 
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-ligV -9r^ + M + 21r^] } z'^ + {-4M + 9r}z^ + 6a^ {5QV + 4r^ - 6r^2 ^ ^ ^_9^2 _^ _^2j } ^ Ar^ 

+ (56g^°r {28r^ + 0^} + 14Q^ {r^[-735M + 298r] + 10[13M - 62r]r'^z'^ + [M- 270r]z^} 

+ Q^{r^ [24955M2 - 19786Mr + 3840r^] - 2r^ [5285M2 - 23710Mr + 8492r^] z'^ + r [763M'^ + 11622Mr 
-4496r2] z^ + 8[-144M + 9O7r]0*'} + {r^ [-26559M^ + 30642M'^r - 11462Mr^ + 1380r^] + [l8645M^ 
-85118MV + 56148Mr^ - 9744r^] z^ - [3097M^ + 5010MV - 4760Mr^ + 1728r^] z'^ + [83M^ + 1886M'^r 
-25348Mr2 + 9792r3] z^ + 2[607M + 198r]z^} - OM^ {[IIM - 4:r]r'^[-2M + r]^ + Sr^ [-AM^ + SOM^r 
-71Mr^ + 16r^] - 9r [4M2 - 19Mr + Sr^] z'^ + [-9M + \^r\z^) + |4M''r [2682r<^ - 2h'i\r^z^ 
+780r2z'* - 55z^] + M^r [8315r^ - 45844r6z2 + 4336r'*z* + 20872r2z^ - 1335z^] - 2M^ [7865r^ - 25324r^z2 
+22h^r^z^ - 740r^z^ + 109z^] - 6M [309r^ - 2719r^z^ + 2639r^z'^ + Alr^z^ - 26z^] + 144r [r^ - Ylr^z^ 
^\%r^z'^ - 42*^] + 84aV^ {VJ^Q'^r + ISM^rp^ + 15Q2 (-ISr^ + z^) + 4rp2] } + {56Q®r [USr^ 
+37^;^] - 7Q^ [3(1525M - 604r)r'^ + 2r2(299M + h%%r)z^ + (-89M + 1692r)z^] + 2Q'^ [3r^ (8491M2 
-6487Mr + \2<^2r^) + r'* (-1918M2 + 16755Mr - 5604r^) + r {-\7hM'^ + 16965Mr - 7578r^) z^ 
+3(-369M + 544r)^''] + 54M2 [r^ (24M2 - 22Mr + 5r^) + r (-4M2 + 45Mr - 20r^) 0^ + (-3M + lOr)^^] 
+Q2 [4m3 (-6812r^ + UObr^z"^ - bUr'^z^ + Alz^) + QM'^r (4892r'^ - 5365r^^;^ - 3902r2^;* + 307/) 
+1080 (r^ - 2r^z'^ - 5r^z^ + 2rz^) - 15M (671r^ - 1839r''2^ + 399r^0'' - 3^;^) p^] } + 6a^ {l68Q*^r [26r^ + \7z^\ 
+56Q'' [r''(-240M + 89r) + 2r'^{-&\M + \3r)z^ + (lOM - 63r)2^] - 45M2 [(13M - 6r)r2 - (M - 8r)z2] p4 
+Q2 [tmV (I459r^ + hh^r^z^ - 45^) - lOM (689r'' - 519r2^^ + 24^) + 360r (2,r^ - 4^^) p^j } j At^ 

+ O(Ar^) 

(5.3c) 



o 02: sm 

ma 



\\2Q^ p^ ^i^^a" + + {20M2 - 28Mr + 9r^} + 2 {16M - 9r} r^^ + 3^ 



6720A2pi8 

+ Ga^ {50^ - lOMr + Gr^ - 4^^} _ jiOMr - Tr^ + Sz^} j At^ 

f (56Qi° {31r2 + 4^2} - 28Q^ {5[88M - 45r]r3 + 2r[4M + \\hr\z^ + 131z^} + {r^ [32515M2 - 32704Mr 

+7912r2] - 2r2 [1757M2 - 18928Mr + 7508r2] z^ + [259M2 + 16128Mr - 9752r2] + 4104z^} 

- 2Q'' {r^ [l8949Af''' - 27967MV + 13136Mr2 - \9^2r'\ + r'"' [-4662M3 + 36765MV - 26924Mr^ + 4896r3] z^ 
+r [581M3 + 10711MV - 14632Mr2 + 4536r3] z'' + 3 [iSTM^ + 2428Mr - 904r2] z*^ - 390z^} 

- 18M2 {[7M - 8r]r3[-2M + rf + 18r2 [GM^ - 7Mr + 2r2] z^ + 3[15M - 8r]rz'* + 2z'^} 

+ Q2 {28M*r2 [599r^ - 210r^z2 + 55z^] + 8MV [-3999r^ + 5881r''z^ + 1135r^z* + 327^^] + [21499r^ 
-47896r*^z^ - 21452r*z^ + 12060r2z'' + 405z*] - 120Mr [SOr'^ - \%\r'^z^ + 62r^z^ + 13z''] p^ + 144 [4r^ - 18r^z^ 
+12r^z^ - z^\ p^} + 84a {l76Q^ + ISM^p^ + 60Q^ [-4Mr + p^] } + 4a^ {28Q^ [l27r^ + 46z^] 
-14Q^ [(1328M - 603r)r^ + 2r(178M + 33r)z^ + 345z*] + 27M2 [2(7M - 5r)(2M - r)r^ + 4(9M - 5r)rz^ 
+5z^] p"* - Q2 [28MV (665r'' + 58r2z2 + 41z'') + ZM"^ (-7755r*5 + 2152r''z2 + 4049r'^z'^ + 190z*^) - 540 (2r^ 
-Sr'^z'' - 2r2z^ + z®) + 30Mr (303r4 - 416r2z2 + 9z'') p^] + [7M'^ (4583r'* + 750r2z2 + 55z^) 
+42Mr (-671r'' + 130r2z2 + 369z'') + 6 (992r'' - 1355r2z2 + 61z*) p^] } + 6a^ {l512Q'^ [3r^ + 2z^] 
-224Q'* [{Q7M - 27r)r^ + 8(5M - 2r)rz2 + llz''] - 90M2 [7Mr - + 5z^] p^ 

[-1120Mr (8/ + ^r^z^ - Zz^) + 7M'^ (l755r^ + 9^Qr'^z^ + 59z'') + 360 (4r^ - ^z^) p^] } ) Ar^ + O(At^) 

(5.3d) 
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2 

q a 



11340g>"A (2a^ + Q - 2Mr + - z^) At 



4 2 

r z 



75600A3/2pi7 

+ (q^ {23807r'' + SGSOSr^z^ _,_ 96322^} - {[119701M - 4900Ar]r^ + 8[21670M - 4053r]r^z^ + r[32923M 

+60726r]z^ + 44146z^} + 2025^/^ {r^[~2M + - 3r[-2M + r]z'^ + 2*} + iQ^ {2M^ [32985r^ + 44155 
+5033r2z'^ + 7752*^] - 6Mr [SBlSr'^ - 3695r'^z'^ - 75492"*] + [9749r'* - 20300r2z2 + 2999z^] p'^} 
- 3Q^ {M^ [37174r^ + 47910r5z2 + 5970r3z* + 4450rz'^] - [40349r^ - 20415r''z2 - 34905r2z'' - 2725z^] p^ 
+9Mr [I509r'* - 3090r2z2 + 2952^] p'^ - 1350 [r^ - 3r^z^ + z^] p^} + 675aV^ {124Q4 + l5M^p^ + 30Q2 [-5Mr 
} + 675aV^ {Q® [I45r2 + 103z2] - 2Q^ [(227M - 8Ar)r^ + {U3M - 30r)rz^ + biz^] - 15M2 [2Mr 



2 + z^] p^ + 2Q2 [M^ {7r^ + z^) (25r2 + llz^) _ 2Mr [b&r'^ - 332^) p^ + 15 (r^ - z^) p^] } j At"* + O(At^) 

(5.3e) 



As expected, in the limit a ^ 0, these reduce to Eqs. (4.3 1 



VI. CONCLUSIONS 

Our previous work in Paper I was only valid provided the particle's path is a geodesic of the background spacetime. 
We have shown here how this may be extended with reasonable ease to allow for non-geodesic motion. Some interesting 
insights have come as a result of this study. 

First, it is interesting to note that, although the covariant approach still has the great benefit of being sufficiently 
general to be applicable to any spacetime of interest, we must convert to a coordinate expansion in order to obtain 
a final result. This is not very suprising as the benefit of using a covariant expansion hinges on having a geodesic 
along which to expand. The coordinate expansion, on the other hand, more naturally characterises motion which is 
not along a geodesic. However, as Sec. |V] demonstrates, beginning with the covariant expansion and computing the 
coordinate expansion from it allows the calculation to be done for space-times which previous methods have not yet 
been able to do. 

It is also interesting to look more closely at the presence of the Christoffel symbols in the relation between the 
covariant and coordinate expansions. In the present context, these take on a particularly clear form by being directly 
related to the 4-acceleration of the particle. In particular, the leading order at which this effect appears is given in 



Eq. (2.11c I, which in this case becomes (since the first two orders, v and Va, are identically zero): 

vttt = vttt + ivtrVlt = ^ttt + 3utra'' (6.1) 

where is the radial component of the 4-acceleration of the particle. It is immediately clear from this expression 
that the deviation from the result for geodesic motion is a direct consequence of the particle's 4-acceleration. 
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